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Abstract 

We propose a new way to determine the squark mass based on the shape of di-jet invariant 
mass distribution of supersymmetry (SUSY) di-jet events at the Large Hadron Colhder 
(LHC). Our algorithm, which is based on event kinematics, requires that the branching 
ratio B{q — )• qZi) is substantial for at least some types of squarks, and that ?ti~ /m? <^ 1. 
We select di-jet events with no isolated leptons, and impose cuts on the total jet transverse 
energy. Elf* = Exiji) + ET{j2), on a = ET{j2)/nijj, and on the azimuthal angle between 
the two jets to reduce SM backgrounds. The shape of the resulting di-jet mass distribution 
depends sensitively on the squark mass, especially if the integrated luminosity is sufficient 
to allow a hard enough cut on Eff^ and yet leave a large enough signal to obtain the mjj 
distribution. We simulate the signal and Standard Model (SM) backgrounds for 100 fb~^ 
integrated luminosity at 14 TeV requiring E^^ > 700 GeV. We show that it should be 
possible to extract m^ to within about 3% at 95% CL — similar to the precision obtained 
using mT2 — from the di-jet mass distribution if mq ~ 650 GeV, or to within ~ 5% if 
mq ~ 1 TeV. 
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1 Introduction 



Supersymmetry (SUSY) is one of the best-motivated extensions of the Standard Model (SM) pQ. 
There have been numerous studies of strategies by which superpartners may be discovered at the 
LHC. The SUSY reach of the LHC is usually expressed in terms of the masses of coloured gluinos 
and squarks, expected to be the most copiously produced sparticles [21 |3l H] . More recently, 
the focus has shifted to how well sparticle properties can be determined at the LHC, with the 
most attention being paid to sparticle masses, the underlying motivation being that knowledge 
of the sparticle spectrum will lead us to the mechanism by which SM superpartners acquire 
their masses. The problem, of course, is that it is not possible to construct mass peaks because 
every SUSY event (in _R-parity conserving SUSY models that we focus upon in this paper) 
includes two undetected neutralinos {Zi), which we take to be the lightest supersymmetric 
particle (LSP). 

The variable Meg = ^^^^^ ExiJi) + constructed from transverse energy (Ex) provides 

a rough measure of the mass of the lighter of gluinos/squarks [5], [3]. The end point of the 
same flavour, opposite-sign dilepton mass spectrum, first mentioned in Ref . [6] , yields a precise 
measurement of the mass difference, m^^ — m^^, between neutralinos. Hinchliffe, Paige and 
collaborators [S] pioneered systematic studies of the information that may be gleaned from 
end-points of a variety of mass distributions in LHC event samples with cuts to select out 
SUSY events over SM backgrounds, and showed that (modulo some discrete ambiguities) in 
favourable cases with two-body decay cascades, it is sometimes possible to also construct the 
masses, rather than just mass differences: in this connection, see also Ref. Since then, 
other techniques have been suggested to obtain sparticle masses. These include the use of the 
mT2 variable [S], its cousins rnxGcn and related variables [1], as well as their generalization to 
asymmetric decays [TU], the so-called matrix element method [TT], the presence of kinematic 
cusps in distributions [12], or through multi-lepton channels [T3] . 

In this paper, we propose a new way of measuring masses for the case of event topologies 
with one step cascade decays of the type: 

p(l) + p(2) ^ ^(1) + /)(1) + ^(2) + Jji2) (1) 

where P^*-* are the parent particles (squarks in our analysis) initially produced in the hard 
scattering, D^'^h are the invisible (lightest SUSY particles, LSP) daughters of P^''\ while each 
q^^^ manifests itself as a jet in the detector. None of the existing mass measurement methods 
allow for a separate determination of both the parent and daughter masses from just the process 
([1]). For instance, the much-studied mT2 method yields a measure of (mp — m^)/2mp which, 
only if m|, / mp ^ 1 , provides a measurement of the parent mass [H] . If more complicated decay 
chains [U] or several production processes [2] are accessible, it may be possible to extract the 
individual parent and daughter masses. The method presented here offers an alternative way 
of determining mp (also in the regime m|,/mp <^ 1) from the dijet invariant mass distribution, 
and relies on completely distinct kinematical features from the standard mT2 procedure. 

In our analysis, we assume that the daughter LSP is a bino-like neutralino Zi, while the 
parent is a right squark, which dominantly decay via — )■ g^'ij^ Thus our signal from right 



If instead, the LSP is a wino-like neutralino, ql would dominantly decay to charged or neutral winos. Since 



squark pair production is exactly two hard jets, no isolated leptons or photons, together with 
missing transverse energy (E'™'''^). It has been pointed out [16], and since confirmed [171 [H] 
that SM backgrounds to this signal can be controlled by requiring acollinear jets, even without 
the use of i?™*^*^. If gluinos are marginally heavier than squarks, they decay to squarks and an 
undetectably soft jet, so that gq production then adds to the squark dijet signal. If gluinos are 
light enough so that squarks dominantly decay to gluinos, we do not have the dijet plus E?^^^^ 
signal, and our analysis does not apply. Aside from our assumption about the mass ordering 
between gluinos and squarks and m^^ ^ nig, we endeavor to leave our analysis as model- 
independent as possible, and at the same time realistic in that we include all complications due 
to contamination from other SUSY processes and SM background. We will see below that in 
addition to yielding rriq our analysis may also serve to approximately constrain rrig. 



2 The acollinear dijet signal 

The cross section for pair production of is determined by SUSY QCD just in terms of rrig and 
rUq^, independent of the details of any model. For the mass ordering of interest to us, the decay 
Qji — )■ qg is kinematically suppressed, so that right squarks dominantly decay to the (bino-like) 
LSP. SUSY contamination to the acollinear dijet signal is, of course, model-dependent. For 
definiteness, we assume gaugino mass unification at the GUT scale, take fi = rrig, tan/3 = 10, 
and use the scalar mass unification condition m? ^ m| -|- 0.7m? as a guide to the slepton mass. 
These choices have little effect on the right squark signal other than restricting the gluino to 
be not much heavier than the squark, since otherwise the slepton becomes too light. 

We use ISA JET v7.78 [TH] for our simulation of the SUSY signal at the LHC. We assume 
a toy calorimeter with cell size Arj x A(j) = 0.05 x 0.05, extending to |?7| = 4. Here ri,(j) 
denote the jet pseudo-rapidity and azimuthal angle. The hadronic calorimetry (HCAL) energy 
resolution is taken to be 100%/^/ET (B 5%, where the two terms are combined in quadrature. 
The electromagnetic calorimetry (ECAL) energy resolution is assumed to be 5%/ ^/ET 0)0.55%. 
We use the cone type ISAJET jet finding algorithm, with a cone size AR = + A02 < 0.4 

to group HCAL energy depositions into jets. Jets are then defined to be hadronic clusters with 
ET^jet) > 50 GeV. Leptons with transverse momentum pt{1) > 5 GeV are defined to be 
isolated if the visible activity within a cone of AR < 0.2 about the lepton direction satisfies 
EE^^"^ < 5 GeV. 

In the analysis of our signal, we require exactly two jets with, 

1. ET{j)> 50 GeV, 

2. azimuthal angle separation A(f)jj < 1.5, 

3. Eif^ = EtUi) + Et{32) > 700 GeV, and 

4. a = ETij2)/mjj > 0.5. 

where ji, j2 are ordered by jet Et (-EtOi) > ET{j2)) and rrijj is the invariant mass of the two 
jets. We veto events with isolated leptons or a third jet with Et^J) > 50 GeV. These cuts and 
veto criteria have been shown to very effectively remove the QCD dijet background [T71 [T8] . 

the SM daughters of the chargino decays would be very soft, the analysis we describe for qfj would then apply 
to Ql- 
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Figure 1: SUSY signal cross section contours for the acollinear dijet signal from supersymmetry 
at a 14 TeV pp collider, after cuts listed in Eq. ([2]). The dashed line shows the boundary of 
sampled area taken approximately to be rrig ~ 0.7mg. The upper left part of the plane is cut 
off because, for our choice of slepton masses, this region is forbidden either because the slepton 
is lighter than Zi or mj drops below 100 GeV. The cross section is largest just above rriq = rrig 
where it can easily exceed 10 fb. 

The SM background is then dominated by {Z — j- uu) + jj events [T71 [18], which has a cross 
section azjj = 6 fb, with an additional contribution of <Twjj = 1-8 fb from W ^ iu events 
where the lepton is not identified. 

Contours of the signal cross section at a 14 TeV pp collider after the cuts listed above are 
shown in the rriq — rrig plane in Fig. [H We cut the contours off at large values of rrig where 
the slepton either becomes lighter than Zi, or when it violates the LEP lower mass limit on its 
mass. The cross section drops off when rrig falls significantly below rriq because squark decays 
to gluinos begin to eat into the branching ratio B{q ^ qZi). The region where we terminate 
the 1 fb contour has recently been excluded at 95% CL by the CMS and ATLAS experiments 
at the LHC [20]. 

When the gluino is heavier than the squark the g ^ q + q decay channel opens up, and gq 
events can pass the Nj = 2 cut if the q forms a very soft jet, below the jet Ex threshold. This 
roughly doubles the signal just above the diagonal rUg = rrtq line. However, this contribution 
falls off with increasing rrig (and fixed rrtq) because the additional jet becomes more readily 
identified in the detector. If instead rrig < rriq, the q —> g + q decay channel takes over quickly 
as the mass difference increases. This decay channel leads to multiple jets and the signal drops 
below 1 fb at rrig ~ f'^g- The signal is also suppressed at low {rriq, rrig) due to the E^^ cut. 
Nevertheless, we see that there is a sizable region where there are several hundred events in the 
SUSY dijet channel for an integrated luminosity of 100 fb~^. 



3 Standard Model backgrounds 

The dominant SM background to the SUSY signal in the acolhnear dijet channel comes from 
the (Z— 7- uu) + 2j events. Because the properties of jets in Z + 2j events at the LHC cannot 
depend on how the Z decays, this background can be directly obtained by scaling LHC data 
on Z ^ e^e~ + fi~^fi~+2 jet events by a factor, 

BR{Z uu) 
^ ~ BR{Z e+e~/fi+fi~) 

The lower rate in the (Z— )■ ii)+jj channel, however, implies that the fluctuations in the inferred 
background will be larger by relative to the corresponding fluctuations in the Z{^ uu) + 2j 
sample. In our analysis we have used the Monte Carlo event generator AlpGen [21] to simulate 
the SM (Z— )■ background, and then scaled the uncertainty in each mjj bin by to 

correctly simulate the corresponding fluctuations when this background is obtained from the 
Z ii + jj data. 

In addition, as already mentioned, W + jj events make a subdominant but sizeable contri- 
bution to the background. Occasionally the charged lepton (most likely a hadronically decaying 
r) from a high px W decay becomes buried in the hadronic jets and/or is undetected because 
it is either too soft or in uninstrumented regions of the detector, while the neutrino is hard 
enough to reduce the angular separation of the two jets. The relatively small W + jj back- 
ground compared to the leading Z + jj one reduces the relevance of -|- jj contamination. 
While the subtraction of the Wjj background is not as simple as for the background from 
Zjj production, we may expect that the W + jj signal can be measured in a control region 
(such as the dijet + 1 lepton channel) and then extrapolated to the signal region (dijet + 
leptons)!! In our analysis, just as for the Zjj background, we assume that we can use AlpGen 
[2T] to simulate and subsequently subtract the Wjj background from the total di-jet sample, 
but include the corresponding statistical uncertainty in the evaluation of when performing 
our fits to extract sparticle masses. 

The reader may well wonder whether it is further possible to enhance the SUSY signal 
relative to background via a requirement on E^^^^. With this in mind, we show the E^^'^'^ 
distribution of SM Z/W + 2 jet backgrounds in Fig. [2], along with the corresponding distribution 
for the two SUSY test cases that we introduce below for our study of how well the squark mass 
may be extracted at the LHC. The total SM background is 7.8 fb after the cuts in Eq. ([2]) 
to be compared with the signal of 16 fb (11 fb) for rriq = 650 GeV (1 TeV). Note that the 
angular and i?^* cuts preclude very low values of E^^'^^. As the signal E^^^^ distribution moves 
to higher values with increasing rriq, it may appear tempting to impose an additional EJ^^^^ cut 
at 750 GeV to further enhance the signal events relative to the background. We have checked, 
however, that the loss in signal statistics largely negates the benefit of a lower background; the 
E'™'^'^ cut does not yield a significant improvement to the determination of rriq. 

^Here, we are assuming that SUSY does not significantly contribute to Zjj events. 

•^This could be complicated by the fact that SUSY events may also make a contribution to the jets plus 
lepton channel, but typically, one would expect that the signal jet multiplicity is larger. 
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Figure 2: The EJp^^^ distribution for the leading SM backgrounds from Z + jj and W + jj pro- 
duction, together with the same distributions for the signal points {mg,mg) = (650,667) GeV 
and (1000, 1030) GeV at a 14 TeV pp collider after the cuts in Eq. ([2]). The Z + jj background 
shown here is the SM prediction for invisible Z decays and does not show the enhanced fluc- 
tuations that would result if this background were extracted from the (Z — e'^e~//i+/i~) + jj 
data as discussed in the text. 

4 Squark mass determination 
4.1 Methodology 

To facilitate the discussion of our strategy for squark mass determination, we show a scatter 
plot of the two jet energies in Fig. [3] for (a) rriq = 650 GeV, and (6) = 1 TeV. Although we 
include all SUSY events in this sample, about 40% (> 90%) of the events come from squark pair 
production in frame (a) where the gluino-squark mass gap, AM, is just 17 GeV, (in frame (6) 
with AM = 30 GeV) with the bulk of the remaining events arising from gqpi pair production. 
In the former case, although squark production seems to be subdominant, we should remember 
that the gluino decays into a squark and a (very soft) jet, so that these contaminating qg events 
have very similar kinematics as squark-pair events. Hence qg events do not interfere with our 
strategy to extract the squark mass and, in fact, help in that they increase the signal sample 
that we can use0 

The gray scale on the scatter plot shows how the events are distributed by invariant mass. 
We see that the bulk of events with low invariant dijet mass nijj have small values of £'7^(72) 
and (because of the E^^ cut) correspondingly large values of -E'r(ji). This is not surprising 
since (except when both jets are very forward) our jet cone algorithm generally requires the 
jets to have a spatial angular separation between them, and precludes the smallest invariant 
masses if ET{j2) is also sizeable. Despite the fact that j2 is soft, these events readily satisfy 

■^If the gluino-squark mass gap becomes substantial, the additional jet fom gluino decay frequently fails the 
third jet veto, suppressing the gluino signal. In case mg ^ rriq the gluino production is suppressed and q qZ\ 
remains the only viable signal. 




Figure 3: A scatter plot of the jet transverse energies for dijet events from supersymmetry 
at a 14 TeV pp collider for Frame (a) ruq = 650 GeV, and Frame (b) rriq = 1 TeV after the 
{a, ^(pjj, -Ey *) cuts described in the text. The gluino mass is chosen slightly higher than the 
squark mass to enhance event rates by including gluino production channels. The invariant 
mass of the dijets is shown by the gray-scale on the dots. Notice that the region with the 
largest Eriji) at the lower part of the plot is populated by events with low invariant mass, 
while events with high values of £'t(j2) and of rrijj preferentially populate the upper part of the 
plot. This leads to the shoulders in the rrijj distribution for the 650 GeV squark case shown 
in the inset frame; this feature, absent for the 1 TeV squark case, would appear if we require 
E*?* > 1 TeV. 



the a requirement because rrijj is also small. Relative to these events with the lowest values 
of nijj, events with somewhat larger dijet masses are possible if either £'T(i2) or the space 
angle between the jets is increased. If nijj becomes larger because £'t(j2) is increased, the 
a cut is still satisfied. However, events with increased values of nijj due to a larger angular 
separation but where ET{j2) remains small have smaller values of a, and so a reduced efficiency 
for passing the a > 0.5 cut. We emphasize that without the Et^^ cut these "intermediate rrtjj 
events" would also be obtained for smaller values of -E'r(ji) and correspondingly larger values 
of -Br(j2) that would readily satisfy a > 0.5. Finally, events with large values of Ej'{j2) close 
to the -Et(ji) = Erp[j2) line naturally tend to have high values of dijet masses (as evidenced by 
the lack of dark points in this region of the scatter plot) as well as of a. Of course, events with 
exactly back-to-back jets are always eliminated by the a cut (and also the A0 cut), which is 
why this cut is an efficient veto against QCD dijet events. 

To recap, these qualitative arguments suggest that, compared to events with very small or 
relatively large values of m^j, events with intermediate values of dijet invariant mass are less 
likely to satisfy a hard cut on Et^^ simultaneously with the requirement a > 0.5. The dip in 
the di-jet invariant mass distribution near rrijj ~ 350 GeV in Fig. ^a) is a manifestation of 
exactly this feature. This deficit of intermediate events, as we will see below, which plays an 
important role in the extraction of rrig from the di-jet data, does not show up in the 1 TeV 
case in frame (b). To better understand this, we have analysed the kinematics of di-jet events 
from squark decay in the Appendix in more detail. In the approximation that the hard jet ji 
is in the direction of its squark parent, we show that with E^^ > 2Eq = (m? — m'~_J/'mq, high 
mass events with large opening angle between the jets, as well as low mass events with a small 
opening angle between the jets readily satisfy that a > 0.5 requirement. We find, however, 
that intermediate mass events where the second jet is roughly perpendicular to the hard jet (in 
the squark CM frame) can never satisfy the a cut if Etf*^ > 2.6Eq. Although our analysis in 
the Appendix neglects events from gq production for which the kinematics is slightly different 
because of the soft jet from gluino decay, and also any QCD radiation, it nevertheless confirms 
the general features that we inferred from our qualitative arguments of this Section. More 
importantly, it fixes the lower limit on E'^* in order for the dip to become apparent in the 
rrLjj distribution to be ~ 2Eq ~ rrig if m^^ ^ nig. We now understand why we do not see a 
corresponding dip in Fig. E] (6) — the -E^* > 700 GeV cut is not hard enough, and the would-be 
dip region is populated by events with not-so-large values of Exiji), and concomitantly larger 
values of £'^(^2) that then pass the a > 0.5 cut. Indeed, we have verified that by hardening the 
E"^* cut to 1 TeV we recover the dip also in this case. We have not shown this because the event 
rate then becomes too small; i.e. despite the appearance of the dip the reduced event rate does 
not lead to any improvement in the determination of rrig, at least for an integrated luminosity 
of 100 fb~^. For this reason, we leave the E^^ cut at 700 GeV throughout this analysis. 

We fit the synthetic rrijj data for the cuts in ([2]) to theoretical templates for a grid of {rUg, rrig) 
in order to extract the squark mass. We use an integrated luminosity of at least 1000 fb~^ to 
generate these templates. We will see below that the squark mass can be extracted with greater 
precision if the integrated luminosity is sufficient for the implementation of a large enough Ei^^ 
cut so that the dip structure in the di-jet mass distribution (which is very sensitive to mg) is 
clearly evident. 



4.2 Results 



We now discuss the precision with which the squark mass can be extracted via fits to the di-jet 
mass distribution. We use the rrig = 650 GeV and rriq = 1000 GeV cases introduced above as 
illustrative examples. In order to extract the squark mass we fit the shape of the signal rrijj 
distribution — obtained after statistical subtraction of the backgrounds from the signal plus 
background sample as explained in Sec. [3] — to the corresponding distributions obtained using 
templates that have been independently generated for a grid of values of {rrig, nig). Given that 
a signal has been detected in the di-jet sample, we define ioi the normalized signal spectrum 

where i sums over all rrijj bins, each one 30 GeV wide, and bins with less than 5 signal 
events are dropped. Here, 0° correspond to the "data" for 100 fb~^ in the bin, while (pi 
is the corresponding expectation for the value of {nig, rrig) from the template obtained with 
an integrated luminosity exceeding 1000 fb~^. The uncertainty 5, = ^JNf"■^ /Nq°^ gives the 
statistical error in the shape of signal spectrum after background subtraction. In this analysis 
we do not include systematics. We minimize over a grid in {rrig, rrig) values, with its 

minimum being the best fit for {nig, nig), and as usual map out contours of constant Ax^ in the 
mq-rrig plane. As mentioned in Sec. [2l gluino events only contribute to the signal if nig ~ rrig. 
The gluino contribution improves the signal statistics but also has a small effect in the rrijj 
shape, and so affects the sparticle mass extraction as will be discussed shortly. 

We note that Eq. (jl]) assumes perfect evaluation of the theoretical spectrum 0j. In practice, 
although we use an order of magnitude larger integrated luminosity for the calculation of the 
expectation from the templates, significant fluctuations remain distorting and even fragmenting 
the la (68% confidence level (CL)) contour. For this reason, we only show results for the 
extraction of sparticle masses at the 2a and 3a contours in the following. 

Fig. m illustrates the results of our fit of the di-jet mass spectrum after the cuts of Eq. ([2]) 
for the two SUSY cases. The inner (outer) line corresponds to the 2a {3a) contour where Ax^ 
= 6.2 (11.8) above the minimum value of x^- The reader may be surprised to see that the 
di-jet distribution also serves to constrain the gluino mass. To understand this, we note that 
starting from the best fit value, as we reduce rrig (keeping rrig fixed), squark decays to gluino 
open up and start to eat into the branching ratio B{q ^ qZi). When the gluino is only just a 
very tiny bit below m^, it decays via g — > qq* , with the virtual squark almost on its mass shell, 
and the additional quark jet being too soft for detection. The virtual almost-on-shell-squark 
then decays via q* — qZi and the fit is almost unaffected because events from squark decays to 
gluinos cannot be kinematically distinguished from events with squarks decaying to the LSP. 
However, as the gluino-squark mass difference increases, (1) gluino decays of squarks become 
more significant and begin eating into the branching fraction for LSP decays, and (2) the quark 
daughter of the squark becomes harder and so more likely to be detected. Indeed, if all q gq 
decays would be vetoed by our jet-multiplicity cut, all that would happen would be a loss of 
statistics, causing the ellipse to widen, and ultimately open up, for lower values of nig. However, 
before this can happen, the three-body gluino decay still has quasi-two-body kinematics, but 



with an off-shell squark lighter than rrig. As a result the template gluino events will be similar 
to squark events with a smaller squark mass; this alters the shape of the corresponding rrijj 
distribution, enabling the ellipse to close at low gluino masses. 

If on the other hand, starting from the best fit value we now go up in rrig (again keeping 
rUq fixed), the gluino decays via g — qq. Initially, the quark from the gluino is very soft, and 
gluinos simply act as a source of on-shell squarks. Thus qg production just leads to better 
statistics, and so lead to a more precise determination of rUq, as we have already noted. As the 
template value of rrig is increased with rrig fixed, the distribution of the daughter squark 
from gluino decay, and hence of its daughter jet, is affected. Moreover, the quark daughter of 
the gluino becomes increasingly more likely to be detected, and events with gluinos are more 
likely to be vetoed. These effects combine and together lead to an increase of Ax^. Note, 
however, that the quark daughter of the gluino will not be significantly harder if the template 
squark mass is also increased along with the template gluino mass, so that the template mass 
gap is still small. In this case, the effect that we have just mentioned is somewhat ameliorated, 
causing the error ellipses to tilt to the right as well as narrow in the squark direction as they 
extend to larger gluino (and squark) massesEI 

We see from Fig. Hlthat this analysis constrains the squark mass at 2a to: 

{mg,mg) = (650,667) GeV : rrig = 635 - 690 GeV , . 

{mg,mg) = (1,1.03) TeV : nig = 935 - 1040 GeV 

The better precision in the first case is partly due to the fact that the higher statistics allows 
us to make use of the dip structure that we discussed in Sec. 14.11 We note also that the gluino 
mass is constrained to lie between 630-760 GeV (940-1220 GeV) in the two cases, with the true 
value being significantly closer to the lower end of the range. 

5 A comparison with squark mass determination using 

In this section we briefly examine how the extraction of the squark mass from the di-jet mass 
distribution compares with the corresponding determination using the mT2 variable [8] that has 
received extensive attention during the last few years. There is no question that, as we have 
already mentioned in Sec. [T]an mT2 analysis can potentially lead to a simultaneous determina- 
tion of gluino and LSP masses [H] if squarks are heavy. Here, we focus our attention on the 
more difficult case that the squark is essentially degenerate with the gluino, in which case the 
kink structure on which the gluino and LSP mass determination is based is less evident [13]. 
For comparison purposes, we focus on acollinear di-jet events with the set of cuts Eq. ([2]) for 
the same test cases in Fig. |H As in Sec. 14.21 we assume that the SM backgrounds (Z + jj and 
W + jj) can be subtracted, but still include these in our analysis of the statistical significance 
of the signal. 

^If however, rrig 3> iriq the rate of gluino events would be less frequent, and squark events for this case would 
differ from those for the test point only by changes in the squark Et distribution and the event rate due to the 
much larger value of rrig. It is conceivable, therefore, that there is another region (well into the grey region in 
the figure) where Ax^ is reduced. We have not investigated whether this docs indeed occur. 




Figure 4: 2o" and 3a contours obtained by fitting the shape of the SUSY signal di-jet mass 
spectrum for (a) nig = 650 GeV, and (6) nig = 1 TeV. As in Fig. [1], the shaded area on the 
upper-left is excluded because the slepton is either lighter than neutralino or below 100 GeV. 
The jaggedness of the contours is due to Monte Carlo statistics. 

The mT2 variable is an extension of the well-known transverse mass (mr) [22] for the case of 
two invisible final state daughters. For the one step cascade decay considered here (see Eq. ([1])), 
we have: 

mT2{^x) = min [max{rrS^\rrSpy\ 

Pt(£'i)+Pt(D2)=E7== 

with 

mf = ^ml + 2(ET(g.)^T(A) - PT(g.).PT(A)) , 

where D = Zi, P = q, is the trial daughter mass and the final state quark masses have 
been neglected. The minimization on the right-hand-side is carried out over the partitions of 
the E™*^*^ vector. Since mj'2{fTix = f^n) < ""^p? the value of m^2^(m£)) determines the parent's 
mass. However, because the LSP mass (m/)) is not known, in general it is not possible to obtain 
the parent mass mp (the squark mass in our case). It can be shown that for any value of the 
trial LSP mass rrix, 

mV^T{^^) = Eo + ^jEl + ml , (6) 

with Eq = (m? — m| )/2mg. Thus for any value of rrix, a determination of m^f^ serves to 
determin^ Eq. Assuming that m| /m| <^ 1, this reduces to a determination of the squark 



^Of course, this simple analysis will be affected by the altered kinematics from gluino events contributing to 
the di-jet sample, and also by QCD radiative corrections. 



mass. In this respect, the information that we get from an analysis is identical to that 

we obtain using the method that we have described in Sec. 14.11 In our analysis we take = 
to obtain, 

m^^''{m^ = 0) = 2Eo ^ nig . (7) 

We show the mT2{jnx = 0) distributions for the SUSY points {mq,mg) = (650,667) GeV 
and (1000, 1030) GeV after the cuts of Eq. ([2]) in Fig. |5l We see that the distributions show 
distinct edges around 660 GeV and 1 TeV for the rriq = 650 GeV and 1 TeV cases, respectively. 
To extract the m^2^ value from the mj'2 distributions we fit a "linear kink function" for a 
selected bin range on either side of the visible edge. The error for the fit includes the signal 
and BG statistical errors added in quadrature, as discussed in Sec J4.2[ With this procedure, 
for the two case studies we obtain at 2a: 

{nig, m-g) = (650,667) GeV : rUg = 640 - 682 GeV , . 

{mg,mg) = (1,1.03) TeV : rUg = 970 - 1030 GeV 

We see that for the lighter squark mass case for which we were able to use the dip structure, 
the error is essentially the same with both methods, whereas for the heavier squark case, the 
mT2 analysis yields a somewhat lower erroi^l. 

The squark mass result for the first SUSY point is shifted towards higher values than the 
true value. This is due to the gq contamination, which tends to mimic qq events with a heavier 
squark mass. On the other hand, the 1 TeV squark case has a much smaller contamination of 
qg events (< 10%) and in this case provides the correct value for rriq. 

Before ending this section, we list the pros and cons of the extraction of the squark mass 
from 771x2 or from the di-jet mass distribution described here. 

• Clearly, the di-jet mass method is inapplicable if the squark is much heavier than the 
gluino, since then the di-jet signal from squark production is very small. In this regime, 
771x2 is clearly superior as it potentially offers the possibility of determining both and 

• If nig ~ nig, and the integrated luminosity of the LHC is large enough to enable us to 
utilize the dip structure in the m^j mass distribution, the two procedures yield the same 
information [i.e. a measure of if m| <C tu'^), and with comparable precision. The 
di-jet mass distribution also offers the possibility of constraining m^. This can then be 
compared with a direct determination of from the independent multi-jet event sample 
using, for instance, 712^3- Consistency of the gluino mass obtained from these independent 
data samples would bolster the case that the new physics is indeed supersymmetry. 

^It has been suggested [23] that if the squark pair is produced with a high transverse momentum (^ niq) 
it may be possible, in principle, to determine both rriq and m^^ separately using the 771^2 procedure. We 
note, however, that under 1% of squark pairs will be produced with these large values of the pair transverse 
momentum, and so do not pursue this any further in our study. 

^We have verified that both the central and error values for rriq is weakly dependent on the choice of the 
trial mass nix- The magnitude of the error is, however, somewhat sensitive to the number of bins on both side 
of the kink that we use in our fit. We have shown a conservative range in Eq.® above. 




Figure 5: The mT2 distribution for (a) {mq,mg) = (650,667) GeV and (b) {mq.mg) = 
(1000, 1030) GeV SUY points in Fig 1 after the cuts in Eq. of the text. The trial LSP 
mass is set to zero. The straight hues show the "hnear kink fit" to the signal, and is used to 
extract m^2^. The dashed histogram shows the distribution for the dominant SM BG (Z 
+ 2 jets). We also show the m^2^ value obtained from the fit and its la error. 

• In our analysis, we have assumed that the squarks that lead to di-jet events (mostly 
un and dji for mSUGRA) are essentially degenerate in mass. If these had a small but 
significant mass splitting, the shape of the rrijj distribution would likely be distorted and 
we would extract some sort of averaged mass. The 7717^2 analysis, on the other hand, relies 
on an end-point (which the mass splitting may smear) and would yield the larger of the 
squark masses. 

• Our method for extracting rrig relies only on the jet energies and angles, whereas 7717^2 
also requires E^^^^. While the jet energy scale is likely to be the largest experimental 
systematic uncertainty for mass extraction from the rrijj distribution, the mT2 procedure 
relies on the entire detector, and so is likely to have different experimental systematics. 



6 Summary 

We have proposed a new way to determine the squark mass using templates to fit the shape of the 
invariant mass distribution of acoUinear dijet events at the LHC. Standard model backgrounds 
are suppressed by cuts on a, /^(t>jj, and veto of isolated leptons. Our analysis is independent 
of many details of the underlying model, and requires only that at least some types of squarks 
have a large branching ratio for direct decays to the LSP. It yields a measure of the squark 
mass if m~ <^ m?. While the signal may have sizable EJ^^^^, imposing a E^^^^ cut does not 
lead to a qualitative improvement on the precision with which the squark mass is determined. 

We emphasize that the squark mass is determined by a fit to only the shape of the signal 
di-jet invariant mass spectrum. If mg ~ rriq, the signal may also get a significant contribution 



from gq production which, hke squark pair production, is fixed by SUSY in terms of just 

rriq and rrig. We do not make use of the magnitude of the signal which would depend on model- 
dependent branching fractions for squark and gluino decay. Indeed, the gluino contamination 
and/or the dependence of the squark Erp spectrum on rrig also allows us to constrain rrig, albeit 
with considerably lower precision than rrig. 

As discussed in Sec. 14.11 and in more detail in the Appendix, the cuts on a and E^^ serve 
a dual role. Aside from background reduction, with appropriate choices a > 0.5 (the natural 
value to discriminate the signal from mismeasured QCD di-jet events) and -E^* > rriq cuts, 
we have shown that signal events with intermediate values of rrijj pass the cuts with reduced 
efficiency compared to events with high and low di-jet masses. The resulting dip structure in 
the di-jet mass distribution (see the inset in Fig. [3]a) is very sensitive to the value of rrig, and so 
increases the precision with which the squark mass can be extracted if the integrated luminosity 
is high enough to allow the implementation of a hard enough cut on Ei^^ . 

We have performed two case studies, each with the gluino just slightly heavier than the 
squark, to investigate the precision with which it is possible to extract the squark mass from 
the dijet distribution. Our results are shown in Fig. |H We see that for the mg = 650 GeV case 
for which the dip structure in the di-jet mass distribution should be visible in LHC experiments 
with an integrated luminosity of 100 fb~^, the squark mass is constrained to lie between 630- 
690 GeV at 2cr, while for a 1 TeV squark for which the 100 fb~^ is not sufficient to utilize the 
dip structure, the corresponding range is 935-1040 GeV. We find that in the former case the 
precision that we obtain is essentially the same as that we would obtain using mT2, but the 
systematics of the two measurements are quite different. 

The method that we have proposed using the analysis of di-jet events from supersymmetry 
to determine the squark mass can be generalized to other beyond the SM theories and used to 
determine the characteristic heavy particle mass. The invariant mass construction relies only 
on kinematical cuts, the assumption of a low mass of the stable (invisible) particle and a large 
enough branching ratio for the two body decay of the heavy particle into a jet and the invisible 
lighter particle. 
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Appendix: The Dijet Invariant Mass 

The di-jet invariant mass in Fig. [3] shows distinct features after the cuts in Eq. Here we 
present a discussion of the kinematics behind these distinct shapes and how they relate to 



To be technically precise, we do assume that iiitra-generational squark mixings are negligible. 



the squark mass. For tractability, we neglect QCD radiation, SUSY contamination and SM 
backgrounds. 

As we will see below, the dip structure in the rrijj distribution that we have discussed at 
length in the text mostly arises from the E^^ = -EtOi) +-E't(J2) and a cuts, hence we shall omit 
the cut to start with. Since rrijj, a and E^^ are invariant under longitudinal boosts along 
the beam axis (the Z-direction) and rotations, we can restrict our analysis to the squark-pair 
center of mass frame, with the momentum of both squarks in the XZ plane: 

gf = 7777,5(1, /3 sin 6', 0, /3 cos 6*) , 

= 7mg(l, -/3 sin ^,0,-/3 cos ^) , (9) 

where /3 is the speed of the squarks, and 7 = l/^/l — The jet energy in the squark rest 
frame is Eq = (m? — m~ )/2mq. Its value in the squark-pair rest frame is changed by the boost, 

but is typically of this order. A hard jet (as required by E^^ > 2Eq) results when the daughter 
quark is roughly collinear with its parent squark. In order to satisfy the hard E^^ cut, we will 
for simplicity, make the approximation that the hardest quark is parallel to its parent squark, 
which fixes its four-momentum to be, 

= 7£'o(l + /3)(l,sin^,0,cose). (10) 

However, Og^g^, the angle between q2 and q2 in the squark center-of-mass frame is not necessarily 
restricted to small values and q2 can be emitted in any direction, although, before cuts, ^^292 — 
still is the most likely value if the squark is significantly boosted. The four-momentum of the 
second quark can be written as: 

^2 = -Eo [7(1 + /3 cos6'o), — 7 sin6'(cos 6^0 + /3) + cos6' sin6'o cos 00, sin 6'o sin0o, 
—7 cos 6'(cos 9o + (3) — sin 9 sin 9q cos 0o] , 

where tt — 9q and 0o are the polar and azimuthal angles of q2 in the q2 rest frame (with axes 
oriented so that q2 is moving along the negative Z-axis), respectively. The q2-q2 angle in the 
squark CM frame {9q^q^) is then related to cos 6*0 through: 

_ (3 + cos 9q 

cos (/52'?2 ~ ^ , Q 1 

1 + p COS 9q 

, while the dijet invariant mass {frtjj) is given by 

rrijj = 727^0(1 + /3)Vl + cos^o • (H) 
We now proceed to investigate the shape of the rrijj distribution under the constraints: 

a > 1/2 and E^^* > 2Eq . (12) 
Toward this end, we consider three distinct cases: 

• A) cos 6^^252 ~ 1 cos6'o ~ 1, corresponding to high rrijj values. 



• B) cosOg^g^ = cos^o = corresponding to intermediate nijj values, and 

• C) cos6'g2q2 ^ —1 cos 6*0 ~ — 1, corresponding to low nijj values. 

Case A). This configuration corresponds to q2 emitted nearly along its parent direction. 
Completely back to back jets {9q = 0) always have a < 1/2 and therefore can never satisfy 
Eq.( !T2|) . Thus, we take = + e, where < e ^ 1. Then a and -E^* become: 



Sin6' /l — /3 COS^ COS0O ^ j r^tot omr /■^"'"/^^- n COS 6* COS i 




a-^-e^Y^ ^— ^ and ^ 2Eo[^/ sin ^ - e (13) 

We see from the expression for a above that the constraint a > 0.5 eliminates events with back- 
to-back (or nearly back-to-back) jets, as does the A0 < 1.5 cut in ([2]) of the text. This is, of 
course, why it so effectively reduces the QCD background [16]. The highest di-jet mass events 
after the cuts nevertheless arise from acollinear hard jets with large opening angle between 
them, though without the a cut this distribution would extend out to even higher rrijj values. 

Case B). This configuration corresponds to the case where q2 is emitted perpendicular to 
its parent direction and consequently also perpendicular to qi. The expressions for a and -E^* 
are: 



a = \j — ^-^ \J cos^ 6 cos^ 00 + sin^ (po and E^* = w ^ ^ -Eo[sin 9+{l—P)\J cos^ 9 cos^ (/)o + sin^ ^q] 



If we impose a > 1/2, we have: 



(14) 



cos^ 6 cos^ 00 + sin^ > , =^ f3 < 1/2 



For 00 = TT or we have: 



E'f = y l^Eo[sme + (1-/3)1 cos^l] 
The above expression is maximum for sin^ ^ = 1/(2 — 2/3 + /3^), which gives: 



for /3 < 1/2. Therefore, for (pQ = n ot {i.e. the quarks, the squarks and the proton beam are 
all in one plane), Eq. f|T2|) can never be satisfied. On the other hand, if we look at extremely 
acoplanar configurations with 0o = ^/2 or 3n/2: 



^t' = \l\^Eo[sme+ (1 - /3)] < y^Eo ^ 2.6^0 



We see that such cases are allowed by our cuts. However, imposing Et^^ > 2Eq requires 
large values of |sin0o| ^ we have verified numerically that, for Case B, E^*^ > 2Eq requires 
I sin 00 1 ^ 0.75 - considerably restricting the phase space for configurations that simultaneously 
satisfy the a and E^^ cuts. This analysis, therefore, confirms the qualitative reasoning in 
Sec. 14.11 of the text for the existence of the dip at intermediate rrijj values. More importantly 
this analysis quantifies the magnitude of the cut on E^^ that is needed for the appearance 
of the dip. Increasing the cut to E'^* > 2.6Eq would completely suppress this configuration 
accentuating the dip even further, but this harder cut would lead to a considerable reduction 
in the statistics. 

Case C). Finally, in this case, q2 is emitted anti-parallel to its parent and consequently 
parallel to qi. This is, of course, an oversimplification, since the two would-be jets would then 
be merged into a single jet by our cone algorithm for defining jets, which is why we write the 
expression for a allowing an angle e between q2 and its parent squark Taking Oq = n — e , with 
< e ^ 1, we find that, 

a ^ ''''^^^"/^ and E*?* ^ 2Eo7 sin 9 , (15) 
(1 + /3)e 

independent of 0o. We readily see that this small rrijj momentum configuration can easily 
satisfy the cuts in Eq. ( fT2l) . We point out that a very large squark boost would make this 
configuration less likely. It is thus possible that a very hard cut on E"^* would eliminate most 
events with very low values of rrijj, causing the dip structure to wash out altogether. 

To see how the semi-quantitative analysis with the three cases that we have just discussed 
stands up to inclusion of all events from squark pair production we show in Fig. [6] a scatter plot 
of rrijj versus cos 6*^252 the squark CM frame, after the a > 1/2 and E^f^ > 800 GeV cuts have 
been applied, but with Oq^q^ and 0o allowed to take all values. Events with an azimuthal angle 
separation A0 larger (smaller) than 1.5 are depicted by stars (dots). The points cluster at high 
rrijj values (cos ^5392 ~ 1), as expected, since the squark boost causes the cos Oq^q^ distribution 
to peak at 9g^g^ = 0. We see that the bulk of these events have A0 > 1.5. On the other 
hand, for intermediate rrijj values {cos Oq^q^ ~ 0), which corresponds to Case B discussed above, 
relatively few events pass the cuts. Once we move to large negative values of cos Oq2q2 (Case C), 
the kinematical configuration once again satisfies the cuts, even though the squark boost tends 
to suppress emission in this direction. We thus see how the double-peaked rrijj profile shown in 
Fig. [3]a arises, once we require E^^ > 700 GeV, a > 0.5 and Acpjj < 1.5 to increase the signal 
cross-section and reduce the SM background. As can be seen from Fig|6l the A0jj < 1.5 cut 
suppresses the high rrijj values but does not change the double-peaked structure. 

Finally, our analysis also makes clear that if we reduce the cut on E^^ and allow events with 
-^T* ~ 2i?o, kinematic configurations of Case B will be allowed for any 0o value and the cut 
suppression of intermediate rrijj values will no longer occur, causing the dip to disappear. We 
then expect the invariant mass distribution to present a single peak at intermediate rrijj. This 
is indeed what happens in Fig|3f> for the {mq,mg) = (1000, 1030) GeV point, where Eq = 490 
GeV and the cuts are > 700 GeV, a > 0.5 and Acpjj < 1.5. 




Figure 6: A scatter plot of the invariant dijet mass {rrijj) versus the quark-squark angle ^^353 
for {mq,mg,'m^^) = (650,667,84) GeV, after the E^^ > 800 GeV and a > 0.5 cuts have been 
applied. The (blue) star points have A(j)jj > 1.5, while the (red) dot points have Acpjj < 1.5 
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